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Abstract 

A Collective coordinate variable for adding a space dependent potential to the sine-Gordon 
model is presented. Interaction of solitons with a delta function potential barrier and also delta 
function potential well is investigated. Most of the characters of interaction are derived analytically. 
We will find that the behaviour of a solitonic solution is like a point particle which is moved under 
the influence of a complicated effective potential. The effective potential is a function of the field 
initial conditions and also parameters of added potential. 
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1 Introduction 

Solitons are important objects in most of the branches of sciences. They are stable against dispersive 
effects, and propagate very similar to classical point-like particles. 

Topological solitons play important role in non-perturbative aspects of quantum field theories. 
They are widely used as models of particles which are generated as nontrivial solutions of nonlinear 
field theories. Skyrmions are solitons which are used as a model of hadrons. Solitons also widely 
appear in models of condensed matter, such as Josephson junction, fluxon, spin chain dynamics 
and so on. Describing active regions of DNA is another example of soliton applications in other 
branches of sciences [U [2] [3] . 

Sine-Gordon equation is a very famous model with well known topological soliton solutions. 
This model widely appears in nonlinear problems. Some of the above mentioned situations are 
examples of sine-Gordon applications. 

Recently, there has been an increasing interest in the scattering of solitons from defects or 
impurities, which generally come from medium properties. The motivations come from both the 
theoretical and the application part of physics. The effects of medium disorders and impurities 
can be added to the equation of motion as perturbative terms [H [5]. These effects also can be 
taken into account by making some parameters of the equation of motion to be function of space 
or time [5J [7] . There still exists another interesting method which is mainly suitable for working 
with topological solitons [HIS]- In this method, one can add such effects to the Lagrangian of the 
system by introducing a suitable nontrivial metric for the back ground space-time, without missing 
the topological boundary conditions . 

Interaction of solitons with defects mainly investigates using numerical analysis and numerical 
simulations. Some analytical models have been presented which are constructed with using suitable 
collective coordinate systems. Analytical models for the sine-Gordon and 4 field theories has been 
presented base on the method of adding the defects through the background space-time metric 
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[TOl HTj . They predict most of the " soliton-potential" behaviours with a very good precision. This 
method can be used for objects that their equation of motion results from a Lorentz invariant action, 
such as sine-Gordon model, 4> 4 theory, CP N model, Skyrme model, Faddeev-Hopf equation, chiral 
quark-soliton model, Gross-Neveu model, nonlinear Klein-Gordon models and so on. 

In this paper a simple collective coordinate system is presented. This system is constructed base 
on the method of adding the defects by making some parameters of Lagrangian to be function of 
space. This method can be used for other nonlinear field theories which their equations of motion 
are not Lorentz invariant as well as Lorentz invariant field theories. 

Therefore a new collective coordinate system for solitons of the sine-Gordon field theory will be 
presented in section 2. Results of this model are discussed and also they will compared with results 
of Ref. [10] in section 3. Some conclusion and remarks will be presented in final section 4. 



2 Collective coordinate variable 

Sine-Gordon model in (1+1) dimensions ^=0,1, is defined by 



L = -0 M ^V-A(aO(l-cos$ (1) 
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where X(x) = Ao + V(x). V(x) is a potential parameter and carries the effects of the external 
potential. Potential V(x) is a localized function which is nonzero only in a certain region of space. 
The equation of motion for the Lagrangian (JTJ is : 

d^<j> + X(x) sin(» = (2) 

This equation has not analytical solution with a general function for V(x). If we take V(x)=0, we 
have usual sine-Gordon equation with the following one soliton solution 



(x,X(t)) = itarT 1 [ cxp ( \/\> X - X — 



(3) 



where X(t) = xq — Xt. xq and X are soliton initial position and its velocity. In the following 
calculations Ao has been set to Ao = 1 . 

Center of the soliton can be considered as a particle if we look at this as a collective coordinate 
variable. The collective coordinate could be related to the potential by using suitable function 
V(x) in equation ([2]). Thus the model is able to give us analytic description for the evolution 
of the soliton center during the soliton-potential interaction. By inserting the solution ([3]) in the 
Lagrangian (UJ with adiabatic approximation [5J [5] we have 

£ = 2 [x 2 - l) sech 2 (x - X(t)) - 2\{x)sech 2 [x - X{t)) (4) 

X(t) remains as a collective coordinate if we integrate @ over the variable x. 

L = J Ldx = AX 2 -8-2 J V{x)sech 2 (x - X(t)) dx (5) 

It is a general equation for any kinds of external potential V(x). If we take the potential V(x) = 
e5(x) §5§ becomes 

L = AX 2 - 8 - 2esech 2 (x - X) (6) 
The equation of motion for the variable X(t) is derived from ^ 

8X - 4esech 2 (X) tanh(X) = (7) 
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The above equation shows that the peak of the soliton moves under the influence of a complicated 
force which is a function of soliton position, soliton velocity and characters of external potential 
V(x). If e > we have a barrier and e < creates a potential well. Equation has an exact 
solution as follows 

X 2 = ^esech 2 (X Q ) + X 2 - ^esech 2 (X) (8) 

where Xq and Xq are soliton initial position and its initial velocity respectively. 
Hamiltonian density is obtainable from Lagrangian (j4|) as follows 

5£ = 2 {x 2 + l) sech 2 - X(t)) + 2\{x)sech 2 (x - X(t)) (9) 

Collective energy density can be calculated with integration of ^ over the variable x. Energy of 
the soliton in the presence of the potential V(x) = eS(x) becomes 

E = AX 2 + 8 + 2esech 2 X (10) 

It is the energy of a particle with the mass of m=8 and velocity X which is moved under the 
influence of external effective potential. Figure 1 presents energy of a static soliton as a function of 
its position under the influence of the potential V(x) = e5(x) with e = 0.5. Because of the extended 
nature of the soliton, the effective potential is not an exact delta function. 
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Figure 1: Energy of a static soliton as a function of its distance from the potential V(x) = e8(x) with 
e = 0.5. 

By substituting X from |8j into (|10[) one can show that the energy is a function of soliton initial 
conditions Xq and Xq only. Therefore the energy of the system is conserved. 
Topological charge is 

Q = [ °°^dx (11) 
J+oo dx 
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where is topological charge density. For the solution ([3]) we have 

^ = 2sech(x - X) (12) 

Collective topological charge density q(X) can be calculated by Integrating (fT2"j) over the variable 
x. A simple calculation shows that it is constant and independent of the soliton position 'X'. 

Some features of soliton-potential dynamics can be investigated using equations © and (fTU|) 
analytically which are discussed in the next section. 



3 Soliton-potential dynamics 

A soliton-barrier system is modeled with e > in Eq. |(5J) or Eq. (fT0|) . Consider a soliton which is 
placed far away from the potential (namely at infinity), with initial velocity Xq. It moves toward 
the barrier and interacts with it. There exist two different kinds of trajectories for the soliton after 
the interaction with the barrier, depend on its initial velocity, which separate by a critical velocity 
u c . In low velocities , soliton reflects back and reaches its initial place with final velocity Uf ~ — Xq. 
A soliton with an initial velocity Xq > u c has enough energy for climbing the barrier, and passing 
over the potential. At the velocities Xq « u c soliton interacts with the potential slowly and spends 
more times near the barrier, but this situation is not a bound state. Presented model predicts all 
of the above features. 

Energy of a soliton in the origin is calculated from ([TO]) as E(X = 0) = 4A 2 + 8 + 2e. Minimum 
value of the energy for a soliton in this situation is E m ; m (X = 0) = 8 + 2e. On the other hand, 
a soliton which comes from the infinity with initial velocity Xq has the energy of E{Xq — oo) = 
4Ag + 8 . It is clear that it can pass though the barrier if E(Xq = oo) > E m i n (X = 0). Therefore 
we have 

u c (X = oo) = (13) 

Now consider a soliton which is placed at initial position Xq < (which is not necessarily infinity) 
with initial velocity Xq. Equation ||8]) shows that the soliton reaches — oo(+oo) with final speed 

A(oo) = \J X^ + ^esech 2 (Xq) if its initial velocity is less (more) than the critical velocity u c (Xq). 

Critical velocity in this situation is not This soliton has the critical initial velocity if its initial 
energy becomes equal to the energy of static soliton at top of the barrier X=0.1n this situation the 
critical velocity is 

u c (Xq) = ^1 (1 - sech*(X Q )) (14) 

All presented models predict same critical velocity for a soliton at the infinity. But critical 
velocity of a soliton at initial position Xq is model dependent. For example following expression 
has been calculated in [TU] 

/ e cosIi 2 Xq — 1 . 
u c (Xq) = J- , 2y ° (15) 
y 2 cosh 2 X + ^ 

Critical velocity of Ref. [10] (Eq. ([T5|) ) is compared with (TT4|) in Figure 2 with e = 0.5. Results of 
both two models are the same at Xq = oo and Xq = 0. But transient behaviour of two models 
between these limits are different. It is clear that we can not fit one of the models on each other 
by choosing an effective potential. 

There exists a return point if soliton initial velocity is less than the critical velocity. This point 
can be found by comparing the soliton energy at its initial position and somewhere that its velocity 
becomes zero. Therefore we have 

1 2 • , 1 



cosh X s top e cosh Xq 
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Figure 2: Critical velocity as a function of initial position predicted in Ref. [TO] (dash line) and our 
presented model Eq. (I14p (solid line). 



Equation (|16p clearly shows a linear relation between cosh 2 x and cosll 2 Xfj ■ There is another 

linear relation between cosh y 1 y - — and ~. Direct simulation and also presented analytic model of 
[10] confirm these relations too. 

A potential well can be constructed with e < 0. Consider a particle moves toward a frictionless 
potential well from infinity. It falls in the well with an increasing velocity and reaches the bottom 
of the well with its maximum speed. After that, it will climb the well with decreasing velocity and 
finally pass through the well. Its final velocity after the interaction is equal to its initial speed. A 
soliton-well system is constructed with changing e into — e in the above equations. But dynamics of 
a soliton-well system is very different from what we have seen in soliton-barrier interaction. There 
is not a critical velocity for a soliton-well system, but we can define an escape velocity. A soliton 
at initial position Xq reaches the infinity with a zero final velocity if its initial velocity is 



Xescape = \/^Sech 2 (X ) (17) 



Predicted value for this situation in Ref. [10] is 



Xescape ~\/2cosh 2 A -f " S) 



Figure 3 presents X escape as a function of soliton initial position using (fT?)) and (|18[) . Predicted 
X escap e at the center of the potential in two models is different. But one can show that we can fit 
this value by defining an effective potential parameter e e ff ec tive, however the transient behaviour 
between the origin and infinity can not fitted, figure 3 also show that predicted escape velocity for 
infinity in two models is zero. 

Consider a potential well with the depth of e and a soliton at the initial position X$ which 
moves toward the well with initial velocity Xq smaller than the X esC ape- The soliton interacts with 
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Figure 3: Escape velocity as a function of initial position predicted in Ref. [10], Eq. ([15]) (dash line) 
and our calculation Eq. fjlTf) (solid line). 



the potential well and reaches a maximum distance X max from the center of the potential with a 
zero velocity and then come back toward the well for another interaction. The soliton oscillates 
around the well with amplitude X max . 



X m ax = sech 1 y sech 2 X - ^-X% (19) 
Maximum distance has been calculated with another model in [TU] as follows 



X max = sech- 1 J sech 2 X - ^X 2 (l - ^sech 2 X ^j (20) 



Figure 4 presents sech 2 X max as a function of Xo using (| 19[) and (|20p with e = 0.5 and Xo = 0.5. They 
are very similar to each other and it is possible to fit one of them on the other in the origin, with 
defining an effective potential parameter e e ff e ctive- 

Trajectory of a soliton during the interaction with the potential, X(t) follows from ((U as 

x{t) dx 

(21) 



*(*=°) ^X 2 + pech 2 X - \sech 2 X 

for soliton-barrier (e > 0) and soliton-well (e < 0) systems. The metric model predicts following 
complicated equation [10] 

f XW dx 
t = / = (22) 

/X(t=0) l( v1 , 2 A cosh 2 X(%+co S h 2 X n ) _ 2 
^0 3 J cosh 2 X (%+cosh 2 X) 3 




Fig lire 4* sGchP 1 Xfjidx cis & function of initial position predicted in Ref. [10], Eq. (|20p (dash line) and 
our calculation Eq. (I19j) (solid line). 

Figure 5 presents soliton trajectory respect to time for a soliton with initial velocity Xq = 0.1 and 
initial position X — —3 in a potential well of e = —0.3, simulated with (|2Tj) and (|22| . Figure 5 
clearly shows very good agreement between presented model and metric model in Ref. [10] in low 
velocities. 

If soliton initial velocity is lower than the escape velocity the soliton oscillates around the well. 
The period of oscillation can be calculated numerically using equation (|21| . Some simulations have 
been done using (f2Tj) and also (|22|) and the results have been compared. The results are different 
as figure 6 presents. One can adjust period of oscillation by using an effective potential parameter 
^effective- But general behaviour of oscillation in two models are not adjustable at all. It is clear 
because presented model comes from different model of adding the potential to the system. But 
comparing of two models show that the results are acceptable with a very good approximation. 

4 Conclusion and Remarks 

An analytical model for scattering of sine-Gordon solitons from delta function potential barriers 
and also potential wells has been presented. Several features of soliton-potential characters were 
calculated using this model. A critical velocity for the soliton during the interaction with a potential 
barrier as a function of its initial conditions and the potential characters has been found. The model 
predicts specific relations between some functions of initial conditions and other functions of final 
state of the soliton after the interaction . An escape velocity has been derived for the soliton-well 
system. Oscillation period of a soliton in a potential well also has been investigated using this 
model. 

Calculated characters have been compared with the results of another analytical model. All of 
the simulations for soliton-barrier and soliton-well systems show the validity of presented analytic 
model. Therefore we can conclude that the presented collective coordinate method is able to explain 
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Figure 5: Soliton trajectory as a function of time for a soliton with initial velocity Xq = 0.1 and initial 
position X = —3 during the interaction with potential well e = —0.3 simulated with presented model 
Eq.dZQ (solid line) and metric model of Ref. [ID] Eq.(|22]) (dash line). 
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Figure 6: Soliton trajectory as a function of time for a soliton with initial velocity Xq = 0.02 and 
initial position X = during the interaction with potential well e = —0.3 simulated with presented 
model (solid line) and metric model of Ref. [TO] (dash line). 
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most of the features of sin-Gordon soliton behaviour during the interaction with a potential. But 
this model (like analytical model [10]) is not able to explain fine structure of the islands of trapping 
in soliton- well system. This phenomenon is a very interesting features of soliton-potential systems. 
It is expected to find an acceptable explanation for this behaviour using a better model with 
suitable collective coordinate method. On the other hand using this method for investigation of 
other nonlinear models in potentials is an interesting subject. 
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